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^ ! Abstract 
m 

lO ■ A general way to construct chain models with certain Lie algebraic 

or quantum Lie algebraic symmetries is presented. These symmetric 
Q\ ', models give rise to series of integrable systems. As an example the 

chain models with A n symmetry and the related Temperley-Lieb al- 
gebraic structures and representations are discussed. It is shown that 
corresponding to these A n symmetric integrable chain models there 
are exactly solvable stationary discrete-time (resp. continuous-time) 
Markov chains whose spectra of the transition matrices (resp. inten- 
sity matrices) are the same as the ones of the corresponding integrable 
5h ' models. 
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1 Introduction 



Integrable chain models have been discussed for many years in statistical and condensed 
matter physics. Some of them have been obtained and investigated using an algebraic 
"Bethe Ansatz method" JTJ, see e.g., for periodic boundary conditions and H for fixed 
boundary conditions. The intrinsic symmetry of these integrable chain models plays an 
essential role in finding complete sets of eigenstates of the systems. 

On the other hand, stochastic models like stochastic reaction-diffusion models, mod- 
els describing coagulation/decoagulation, birth/death processes, pair-creation and pair- 
annihilation of molecules on a chain, have attracted considerable interest due to their 
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importance in many physical, chemical and biological processes ||. Some of these stochas- 
tic models can be "exactly solved", see e.g., f|. The theoretical description of stochastic 
reaction-diffusion systems is given by the "master equation" which describes the time 
evolution of the probability distribution function ||. This equation has the form of a 
heat equation with potential (i.e., a Schrodinger equation with "imaginary time"). If an 
integrable system can be transformed into a stochastic reaction-diffusion system, e.g., by a 
unitary transformation between their respective Hamiltonians, looked upon as self-adjoint 
operators acting in the respective Hilbert spaces, then the stochastic model so obtained 
is exactly soluble with the same energy spectrum as the one of the integrable system [0. 

In this paper, by using the coproduct properties of bi-algebras, we present a general 
procedure for the construction of chain models having a certain Lie algebra or quantum 
Lie algebra symmetry with nearest or non-nearest neighbours interactions. The models 
obtained in this way can be reduced to integrable ones via a detailed representation of 
the symmetry algebras involved. As an example we discuss integrable models with A n 
symmetry. These models turn out to have an additional Temperley-Lieb (TL) algebraic 
structure, in the sense that the Hamiltonians give rise to unitary representations of the 
TL algebra and can be expressed by the elements of the TL algebra. What is more, we 
find all these models can be transformed into both stationary discrete-time and station- 
ary continuous-time Markov chains (discrete reaction-diffusion models, see e.g.[[H), whose 
spectra of the transition matrices resp. intensity matrices are the same as the ones of 
these A n invariant integrable models. 

In section 2 we first recall the basic properties of bialgebras and then describe their 
use in the construction of chain models with a Lie algebra symmetry resp. a quantum 
Lie algebra symmetry. In section 3 we discuss the A n symmetric integrable models in 
the fundamental representation of A n . We give the representation of a TL algebra in 
these models. In section 4 we prove that these A n symmetric integrable chain models 
can be transformed into both continuous-time and discrete-time Markov chains. Some 
conclusions and remarks are given in section 5. 

2 Chain Models with Algebraic Symmetry 
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2.1 Bi-algebra 

Let A be an associative algebra. A is said to be a bi-algebra if it contains two linear 
operators, the multiplication m and the coproduct A. 

The operation of multiplication m is defined by: 



m: A <g> A -> A 

m(a <S> b) = ab, Va, 6 G A. 
Let {ej} be the set of base elements of the algebra A. Then ([[]) means 



(1) 



m(e« <g> e^) = m£-e fc . (2) 

k 

The tensor mf- determines the properties of m completely. 
The multiplication m is associative, 

m(m ® id) = m(id <8> wi), (3) 

i.e., 

it it 
where id denotes the identity transformation, 

id : A —> A, id(a) = a. 

However in general m is not commutative, i.e., we have 

mop^m, (5) 

or equivalently, m|j- 7^ m* i5 where p is the transposition operator, 

p: A® A® A, p(a®b) = (b®a), Va, b 6 A (6) 

The coproduct operator A maps ^4 into A £g> A: 

A : A -> A® A, 

( 7 ) 



The properties of the coproduct are described by the tensor \i\ . 
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A is an algebraic homomorphism, 

A(ab) = A(a)A(6), Va,6 e A, 

E /4V, = E ™fcp<A*Vf , (8) 

where A (a) and A (6) belong to A <g> A and the multiplication of tensors is defined by 

(ai <g> a 2 )(&i <g> b 2 ) = ai&i <8> a 2 b 2 , Oi, a 2 , 61, fe 2 £ A 
The coproduct is associative 

(A® id) A = (id® A) A, 

(9) 



"1 ) 



but in general not co- commutative 

poA^A, / (10) 

The operation A preserves all the algebraic relations of the algebra A. It gives a way to 
find representations of the algebra A in the direct product of spaces. 

If a bi-algebra has in addition unit, counit and antipode operators, then it is called a 
Hopf algebra. Lie algebras are Hopf algebras with A co-commutative. Quantum algebras 
are Hopf algebras that are not co-commutative, see e.g. |§ and references therein. 

2.2 Chain Models with Lie algebraic Symmetry 

Let A be a Lie algebra with basis e = {e a }, a = 1, 2, n, satisfying the Lie commutation 
relations 

[e a ,e p ] = C2 p e y , (11) 
where C^p are the structure constants with respect to the base e. 

Let A (resp. C(e)) be the coproduct operator (resp. Casimir operator) of the algebra 
A. We have 

[C(e),e a ]=0, a = 1,2,. (12) 
The coproduct operator action on the Lie algebra elements is given by 

Ae a = e a <g> 1 + 1 <g> e a , (13) 



1 stands for the identity operator. It is easy to check that 

[Ae Q ,A e/3 ] = C^Ae T 
From the properties of the coproduct, AC(e) is a rank two tensor satisfying 

[AC(e), Ae Q ] = 0, a = l,2,...,n. (14) 

Let IF denote an entire function defined on the (L + l)-th tensor space A ® A (g> ... ® A 
of the algebra A. From fllil) we have 

[iF(AC(e)),Ae a ] = 0, « = l,2,..,n. (15) 

We consider "a chain with L + 1 sites", i.e., the set {1,2, ...,L + 1}. We call it 
an algebraic chain in the following. To each point i of the chain we associate a (finite 
dimensional complex) Hilbert space Hi. We can then associate to the whole chain the 
tensor product Hi ® H 2 <8> ... (g> Hl+i- For simplicity we use subindices i,j,k... for the 
points in the chain sites. 

The generators of the algebra A acting on the Hilbert space H\ ® i?2 <£> ■■• ® 
associated with the above chain are given by 

£ Q = A L e Q , « = l,2,...,n, (16) 

where we have defined 

A m = ( id ® . . . (g) id ®A)...(id id A) (id A)A, VmelN. (17) 

m times 

E a also generates the Lie algebra A, 

[E a ,Ep] = Cl E 7 . 

We call 

H = Y / F(AC(e)) i>i+1 (18) 

i=l 

the (quantum mechanics) Hamiltonian associated with the chain {1,2,...,L + 1} and given 
by real entire function IF. Here iF(AC(e))j i j + i means that the rank- two tensor element 
lF(AC(e)) is on the sites i and i + 1 of the chain, i.e., 

F(AC(e)) iii+1 = li ® ... ® ® F(AC(e)) ® l l+2 <g> ... ® l L+ i. 
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[Theorem 1]. The Hamiltonian H is a self-adjoint operator acting in H\®H2®---®Hl+i 
and is invariant under the algebra A. 

[Proof]. That H is self-adjoint is immediate from the definition. To prove the invari- 
ance of H it suffices to prove [H, E a ] = 0, a = 1, 2, n. From formula ( |13D E a in (|16D is 
simply 



E a = E( e ") 



(19) 



i=l 



where 



Obviously 



1, 



i+1 



kL+1- 



[iF(AC(e)) iii+1 , (e a ),-] = 0, Vj ^ i, i + 1. 
By using formula (JJUD and (|l"3|) we have 



(20) 



E^( ACr ( c ))wi.E(Ca)i 
i=l j=l 



i-l L 

e a)j + ( e a)i + ( e a)i+l 

j=l j=i+2 



E 

i=l 
L 

e a)i+l] 

t=l 
L 

E[^(AC(e)) M+1 ,(Ae Q ) 



i=l 



By formula (|T5|) we get [if, E^] = 0, a = 1, 2, n. ■ 

The Hamiltonian system given by (|18| ) describes nearest neighbours interactions. Sys- 
tems with A-invariant Hamiltonians and non-nearest neighbours interactions can be con- 
structed by iterating the application of the coproduct operator to the Casimir operator. 
For instance, the following Hamiltonian describes a system with N + 1(N < L) sites 
interactions: 

L-N+l 

(21) 



L-N+l 

H N = E IF{A N C{e)\ t+l _ i+N) 

i=l 



with A^ as in definition (17). H N commutes with the generators E Q , a = 1,2, ...,n, of 



rN 
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the algebra A on the Hilbert space Hi <g) ... <g) since 
"l-jv+i l 



E iF(A^C(e)) M+1) ... )i+iV , ]T(e Q ). 

i=l i=\ 
L-N+l I" i-1 i+AT L 

i=l j=l j=i j=i+N+l 

L-N+l I" i+Af 

£ iF(A Ar C(e)) M+1 ,...,^,E(e Q ), 

i=i [ j=j 

L-iV+1 



i=l 



where the relation A N ([lF(C(e)),e a \) = [A N (lF(C(e))), A N (e a )} = has been used. 



2.3 Chain Models with Quantum Lie Algebraic Symmetry 

Let e = {e a , f a , h a }, a = 1, 2, n, be the Chevalley basis of a Lie algebra A with rank 
n. Let e' = {e' a , f a , h' a }, a = l,2,...,n, be the corresponding elements of the quantum 
(q-deformed) Lie algebra A q . We denote by r a the simple roots of the Lie algebra A. The 
Cartan matrix (a a p) is then 



da/3 = j-( r «' r /3)> d a =^(r Q -r a ). 



(22) 



We introduce a complex quantum parameter q, such that q da ^ ±1,0. The quantum 
algebra generated by {e' a , f' a , h' a } is defined by the following relations: 



Kih'f,] = o, 



IK, fs\ 



K, f'p] 



—a. 



(23) 



together with the quantum Serre relations 



l-a a g 

E (-i) 7 

7=0 
1— V 

E (- 1 ) 7 

7=0 



1 — a a/3 
7 

1 — a af3 
7 



(24) 



where for m > n G IV, 



m 
n 

[n\ q 



\m\ 



[n] q \[m - n] q V 
[n) g [n-l) g ...[2) q [l) g , 



q n — q n 



q-q 



-i 



The coproduct operator A' of the quantum algebra A q is given by 



A'h' c 
AV r 
A7, 



h'®l + l®ti„ 



a 



g -d a h' a + q d a h' a ^ 



/; <8> g" d ^ + q daK ® f a . 



(25) 
(26) 
(27) 



It is straightforward to check that A' preserves all the algebraic relations in (|23|) and ( P4[) . 
Let C q (e') be the Casimir operator of A ? , i.e., [C g (e'), a] =0, Va G A r For any entire 



function IF of C q (e'), we have 



[F(C ff (e')),a] = 0, Vaei, 



and 



[A'F(C 9 (e')), A'a] = 0, Va G A ? . 
Especially, by formula (p5|) one gets 

„±d a h' a _ ±d a h' a ±d a h' a 



Hence 



[A'iF(C 9 (e')),AV 



[A'iF(C g (e')),g 



(2* 



(29) 



(30) 



(31) 
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The generators of A q on a chain with (L + l)-sites are given by 



L+l 



EL 



i=i 



A ,L e 



L+l 



l 

ai 



E (q dah ' a )i ® - ® (q dah ' a )i-i ® (Oi ® ® ... ® (32) 



L+l 



E(^)i ® - ® ® ® (<r dc ^)m ® ... ® fa-*"*-) 



L+l- 



i=l 



[Theorem 2]. The chain model defined by the following Hamiltonian acting in H\ 
g> i?L+i is invariant under the quantum algebra A q : 



i,i+l- 



(33) 



i=l 



[Proof]. Using @ and we get 



L 

E 

i=i 

L 



L+l 



(A'iF(C 9 ( e '))) M+ i,E(^ 
i=i 



E [(A'^co) Wi, TO* + TOwl 



i=i 

L 



J2[(A'F(C q (e'))),A'(h' a )}. i+1 = Q. 



i=l 



From formulae (52) and (3~I|) we have 

[H g ,E' a ] = [E(A'iF(C g (e'))) M +i, 
.i=i 

L+l 

J2(q dah ' a )i <8> .- ® (/ a ' l ") i _i ® (e'J,- ® (q^Sj+i 
j=i 

= £ [(A'JF(C 9 (e'))) M+ i, &)< ® (9-^)^.1 + (g^*-) 
i=i 

Using formulae (EB) and (53) we get 

= E [(A'JF(C g (e'))),A'(e' a )]. i+1 = 0. 



H+l 



I L+l 



i=l 
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Similarly we have 

[H q ,F> a ] = \j:(A'F(C q (e'))) hl+u 
a=i 

L+l 

Y,(Q dah ' a h ® - ® ® ® Gr^W ® ...(g-***-)^ 

= E [(A'iP(C g (e'))) M+ i, (O* ® (<T***4m + C^)* ® 

i=l 
i=l 

Therefore H q commutes with the generators of A q for the chain. ■ 

The Hamiltonian ( p3|) stands for a system with nearest neighbours interactions. Gen- 
erally by using the coproduct operator A' we can construct models with N + 1(N < L) 
sites interactions: 

L-N+l 

(34) 



L-N+l 

H?= F(A' N C q (e%, i+1 _ i+N , 
i=i 



Taking into account the relation 



A rN ([lF(C q (e')),q ±d ° h °}) = [A' N F(C q (e')), q 



,/\\ „±d a h' 



q ±d a h' a i = Q ; 



N times 



we can prove that the Hamiltonian H q has the symmetry of the algebra A q . In fact: 

[H q , H' a ] 



L-N+l 

E 

i=l 



L+l 



L-N+l 

E 

i=i 

L-N+l 



(A'^(C 9 (e'))) Vi+ i,...^,E(^),- 

i+N 

,i+l,...,i+N> 

j=i 



E [(A' N lF(C q (e'))),A' N (h' a 



i,i+l,...,i+N 



L-N+l 

,i+l,...,i+N, 

i=l 
i+N 



Y,{l daK )i ® - ® ((Z^i-i ® (e'Ji ® (g~" Q H+i ® ® {q- daK )i+N 
j=i 

L-N+l 

E [(A'"*(C fl (0)),A'Vo 



i=l 



i,i+l,...,i+JV 
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E (A' iV iP(C g ( e '))) i ,i+l,...,i+JV> 
. »=1 

i+JV 



E(^)* ® ... ® (g da H-i ® (/a)i ® (<r da H+i ® - ® (q~ daK )i +N 

j=i 

L-N+l 

E [(A'"F(C ff («0)),A'V« 



i=i 



i+l....,i+JV 



0. 



The Hamiltonian system ( J33|) is expressed by the quantum algebraic generators e' = 
(h' a ,e' a , fa)- Assume now that e — > e'(e) is an algebraic map from A to A q (we remark 
that for rank one algebras, both classical and quantum algebraic maps can be discussed 
in terms of the two dimensional manifolds related to the algebras, see ||10|| ). We then have 



L-N+l 

H?= E F^' N C q {e\e))\ l+1 _ l+N . 

i=l 



(35) 



In this way we obtain Hamiltonian systems having quantum algebraic symmetry but 
expressed in terms of the usual Lie algebraic generators {e a }. 



3 Integrable Models with A n Symmetry 
3.1 Quantum Yang-Baxter Equation 



The quantum Yang-Baxter equation (QYBE) [11] is the master equation for integrable 
models in statistical mechanics. It plays an important role in a variety of problems in 
theoretical physics such as exactly soluble models (like the six and eight vertex models) in 
statistical mechanics |T2]] , integrable model field theories [I3|, exact S-matrix theoretical 



models f!4|| , two dimensional field theories involving fields with intermediate statistics 
conformal field theory and quantum groups 0. In this section we will investigate 
the integrability of the chain models having a certain algebraic symmetry constructed in 
section 2. We also present a series of solutions of the QYBE from the construction of 
integrable models. 

Let V be a complex vector space and R be the solution of QYBE without spectral 
parameters, see e.g. 0. Then R takes values in End<${y ® V). The QYBE is 

R12R13R23 = -R23-Rl3-Rl2- (36) 
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Here Rij denotes the matrix on the complex vector space V <S> V <S> V, acting as R on the 
i-ih and the j-th components and as the identity on the other components. 

Let R = Rp, p as in (§). Then the QYBE (^) becomes 

-R12-R23-R12 = R23R12R23, (37) 

where R\ 2 — R ® 1, -R23 — 1 ® R and 1 is the identity operator on V. 

A chain model with nearest neighbours interactions having a (quantum mechanical) 
Hamiltonian of the form 

L 

H = £(WW ( 38 ) 

i=l 

is said to be integrable if the rank- two tensor operator 7i satisfies the QYBE relation ([37]), 
i.e., 

(H) 12 (H) 23 (H) 12 = (H) 23 (H) 12 (H) 23 , (39) 

where 

{H) 12 = H®1, (H) 23 = 1®H. 

The Hamiltonian system (^) satisfying relation ( |39"D can in principle be exactly solved 
by the algebraic Bethe Ansatz method, see e.g. JT|. 

3.2 Integrable A n Symmetric Chain Models 

The integrability of the models having a certain algebraic symmetry presented in section 
2 depends on the detailed representation of the corresponding symmetry algebra. In 
the following we investigate the integrability of chain models with nearest neighbours 
interactions and Lie algebraic symmetry A n . 

Let (ciap) be the Cartan matrix of the A n algebra. In the Chevalley basis the algebra 
A n is spanned by the generators {h a , e a , f a }, a = 1, 2, n, with the following algebraic 
relations: 

[h a , hp] = 0, 
[h a , e« = a a p€R, 

(40) 

[ha, fp\ — —Uapfp, 
[ e a, fp] = fiapha, 
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together with the generators with respect to non simple roots, 

e a .../3 7 = [e a , [e/3, e 7 ]...j, f a ...py = [/«, [/a, / 7 ]-..j. (41) 

Let i£ Q/ g be an (n + 1) x (n + 1) matrix such that (E a p) r s = 5 a ~,5/3s, i.e., the only non 
zero element of the matrix E a p is 1 at row a and column j3. Hence 

EapE-fS = Sp^EaS (42) 

and 

[E a p, E y s] = S/3^E a s — 5s a E/3~/. 

For the fundamental representation we take the basis of the algebra A n as 

h a = E aa — E a+ i ia+ x, a = 1, 2, n 

e = {E aP } ] (43) 
> (3 > a— 1,2, ...,n 

/ = {JW J 

Both {e Q } and {/<*} have a total of n(n + l)/2 generators. 



With respect to the basis fl43|), the Casimir operator of the algebra A n is given by 

n(n+l)/2 n 

C An = (n+1) 51 (e a / Q + /«e a ) + ^ a(n + 1 - a)/i^ 



n ra— a 



(44) 



+ X 2a ( n + 1 - " - (3)h a h a+p - a, 

a=l ,3=1 

where a is an arbitrary real constant. 

The coproduct operator A is given by 
A(l) = 1® 1 

A(/i a ) = h a <8> 1 + 1 <§) fo Q , a = 1, 2, n 

A(e /3 ) = e p ® 1 + 1 <g> 1 

/3 = l,2,...,n(n + l)/2, 

A(//j) = //J ® 1 + 1 ® /„ J 

where the identity operator 1 is the (n + 1) x (n + 1) identity matrix. 



(45) 



By flU) and (g5|) we have 



AC An = C An ® 1 + 1 <g> C A „ - al ® 1 

n(n+l)/2 n 

+(n + l) 51 (e a ®/ Q + / a ®e a ) + ^a(n + l-a)/i a (g)/i Q 

a=l ct=l 

?i ri— ct 

+ X X "( n + 1 - a - /3)(^a ® fra+/3 + ^a+/3 ® h a) ■ 
a=l (3=1 

14 



(46) 



It is easy to check that under the representation ( K5| ) Ca„ is equal to n(n + 2)1. 
Therefore the sum of the first two terms on the right hand side of fl46f) is 2n(n + 2)1 x 1. 
In the following we take a in ( fHq) to be 2n(n + 2) so that the terms that are proportional 
to (n + l) 2 x + l) 2 identity matrix will disappear in fl4"6|) . 

From ([|3|) and (|46|) we have 



AC 



n+1 

(n+1) ^ E a p®E Pa 



n 



a(n + 1 - a)(£ Qa - E a+ x >a+ x) ® (E aa - E a+1>a +i) 



a=l 
n n—a 



+ a ( n + 1 — « — |3)[(£oa — E a+1>a+1 ) (g (E a+l3ta+/3 — E a+ p + i^ a+ p + i) 

a=l (3=1 

+ {E 

a +/m ) ® (£ aa -£ a+1 , a+1 )]. 

(47) 

AC J 4 n in ([|7]) is an (n + l) 2 x (n + 1) 2 matrix. Its matrix representation is 
(AC^J^ = S a p[(n + l)5 Qi /( n+ i)+/+i - 1] 

+ (n+ l)[^aj( n +2)+fc+2^/3,(i+l)(n+2)+fc(n+l) (48) 



+<5 



pi, 



f3,j(n+2)+k+2Va,(j+l)(n+2)+k(n+l)], 



where a, (3 = 1, 2, (n + l) 2 , Z = 0, 1, n, j = 0, 1, n — 1, fc = 0, 1, n — j — 1, 
^a,i(n+2)+fc+2 = if a ^ j (n + 2) + k + 2 for all possible values of j and fe. Here we give 
explicitly, as examples, the matrix representations of AC An for n — 1, 2, 3: 



AC Al 



/ 1 • • ■ \ 

• -1 2 • 

■ 2 -1 • 

\ ■ • • 1 j 



(49) 
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/ 2 • 



ACV 



-1 • 3 • 
. _1 . . 

3 • -1 • 

2 



-1 • 3 
• -1 • 
3 • -1 



2/ 



(50) 



/ 3 • 



ACV 



-1 • • 4 
. _i . 

-1 • 

-1 



-1 



-1 • 



V 

where for simplicity • stands for 0. 

[Lemma 1]. ACa„ satisfies the following relation 



-1 
• -1 

. _i . 

4 • • -1 



(51) 



(AC A J 2 + 2AC An - n(n + 2)1 ® 1 = 0. 



(52) 
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[Proof]. From ( fjg|) we have 

(n+l) 2 
(n+l) 2 

= X! [<W(( n + l)<W(n+l)+m — 1) + (n + l)(5 Q ,j( n+ 2)+fc+2<5' / a,(j+l)(n+2)+fc(n+l) 
0=1 

+^/3,i(n+2)+fc+2^a,(i+l)(n+2)+fc(n+l))]- 

[(5/3 7 ((n + l)^^/( n+ i) + i/ + i — 1) + (n + l)(5/3J'(n+2)+fc'+2^7,(j'+l)(n+2)+A:'(n+l) 
+^ 7 ,i'(n+2)+fc'+2^,(j'+l)(n+2)+fe'(n+l))] 
= 5cry[(^ + l) 2 <5a,Z(n+l)+Z+1^7,Z'(n+l)+Z' + l] 

— (n + l)(^ a ,i(n+l)+«+l + &y,2'(n+l)+2'+l) + 1] 

-2(n + l)[<5 ,j'(n.+2)+fc+2^7,(j+l)(n+2)+A:(?i+l) + ^7,j(n+2)+A:+2^a,(j+l)(n+2)+A:(n+l)] 
+(n + l) 2 [5 a j( n+ 2)+k+2^-y,j{n+2)+k+2 + 5a,(j+l)(n+2)+fc(n+l)^7,(i+l)(n+2)+fe(n+l)] 

= 8 aJ [(n — l) 2 5 a ,i( n +i)+/+i] 

— 2(n + l)[5 ,j( n+ 2)+fc+2^7,(j+l)(n+2)+A;(?i+l) + ^7j(n+2)+fc+2^a,(j+l)(n+2)+fc(n+l)] 
+ (n + l) 2 5 a7 [5 a j(„ + 2)+fc+2 + <^Q,(i+l)(n+2)+fc(n+l)] 

a, i(n+2)+fc+2 + 5a,(i+l)(n+2)+fc(n+l)] 

+(n + l) 2 5c, 7 5 ai /( n+ i) + ; + i — <5 a7 
= -2(AC A J a7 + (n + 1) 2 <L 7 - 5 ai 
= -2(AC An ) ay + n(n + 2)5 ai , 
where the identity 

"a,£(n+l)+2+l a,i(n+2)+fc+2 a, (j+l)(n+2)+fc(n+l) — 1, (53) 

I = 0,1, n, j = 0, 1, ...,n- 1, & = 0, 1, ...,n — j — 1, has been used. ■ 
[Lemma 2]. The coproduct of the A n Casimir operator AC Ari has the following prop- 

(AC An ®l)(l®AC An )(AC An ®l) 

-n[(l <g> AC A J(AC An ® 1) + (AC An ® 1)(1 <g> AC A J] (54) 
+ (n 2 - l)(AC An ® 1) + n 2 (l ® AC^J + n{\ - n 2 )l ® 1 ® 1 = 

and 

(1®AC A J(AC A „®1)(1®AC A J 

-n[(l ® AC A J(AC An (g) 1) + (A04„ ® 1)(1 ® ACaJ] (55) 
+ (n 2 - 1)(1 ® AC A J + n 2 (A^„ ® 1) + n(l - n 2 )l ® 1 ® 1 = 0. 
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[Proof]. By using the representation of ACa„ in (^3) we have 

(AC An ® l)aj3 = (ACA n )(a- 7 )/(n+l)+l,(/3- 7 )/(n+l)+l 
= bap[{n + l)5 a - 7 ,Z(n+l)(n+2) ~ 1] 

+ (n + l)[5 Q ,_ 7i ( n+ i)(j( n+ 2) + fc + l))5/3-- 7 ,(n+l)(j(ri+2)+(A:+l)(n+l)) 



(56) 



+5 



/3-7,(n+l)0'(n+2)+fc+l)^a-7,(n+l)(i(n+2)+(fc+l)(n+l))j 



and 



:i ® Ac 



(AC , A n ) Q --(„+l)2(y_i) ij a_( n+ i)2(y_i) 

= £a/3[(n+ l)^Q-(n+l) 2 (7'-l),«(n+l)+'+l) — 1] 

+ (n + 1) [5 Q -(„+l)2(7/_i) ): ,-(„ + 2)+fc+2))^/3--(n+l) 2 (7'-l),(j+l)(n+2)+fc(n+l) 
+^/3-(n+l) 2 (7'-l)J(n+2)+fc+2^Q-(n+l) 2 (7'-l),(i+l)(n+2)+A;(n+l)]) 

(57) 

where a, (3 = 1, (n + l) 3 , Z = 0, 1, n, j = 0, 1, n — 1 and k — 0, 1, n — j — 1 as 
in formula (pEBD, 7 = 1, n + 1 such that (a — j)/(n + 1) and (/3 — j)/(n+ 1) in ( p^ ) are 
integers and 7' = 1, n + 1 in (|57j). 

Using the formulae floljD and ( |57D one can get ([53]) and Q5Ti| ) from straightforward 
calculations. ■ 



From Theorem 1 we know that the following Hamiltonian is invariant under A n 

H = f>(ACUJ M+1 . 



(58) 



For the given representation ( [43"D of /l n the integrability of (f48[) depends on the form of 
the entire function IF. Due to the relation fl52|) in Lemma 1, (ACa„) 1 , I > 2, can be 
expressed as cACA n + c'l ® 1 for some real constants c and c'. Therefore -ZF(ACa„) is a 
polynomial in AC^ up to powers of order two. 

[Theorem 3]. The following A n invariant Hamiltonian is integrable 



n(n+l)/2 n 
n + ({ e a)i{fa)i+i + (/ a )i(e«)i+i) + a ( n + 1 ~ a )(K)i{K)i+l 



1=1 
L 

E 

i=l 



a=l 



a=l 



n n— a 



+ a ( n + 1 - « - /5)((^a)<(^a+/9)<+l + (^a+^)i(^a)i+l) 

a=l /3=1 



(59) 



where Ti = AC An + 1 and the number 1 should be understood as the identity operator, 
on the tensor space Hi eg) ... ® Hl+i (L + 1 is the number of lattice sites of the 

chain) . 

[Proof]. What we have to prove is that TC satisfies the QYBE (|39p, i.e., 

(H) 12(^)23(^)12 = (W) 23(^)12(^)23; 

where 

(W)i 2 = (AC in + 1) ® 1, (W) 23 = 1 ® (AC An + 1). (60) 
From (|60|) we have 

12 

= (1 ® AC An + AC An 1+ (AC An ® 1)(1 ® AC A J + 1 ® 1 ® l)(W)i2 
= (1 <g> AC A J(AC An <g) 1) + (AC An ® l)(AC An <g> 1) 

+(AC An ® 1)(1 ® ACUJ(ACU„ ® 1) + (AC An ® 1)(1 ® ACaJ 

+2AC An Ol + l® AC An + 1 ® 1 ® 1 

and 

(W) 23 (W) 12 (W) 23 

= (W) 23 (l ® AC An + AC An ® 1 + (AC An ® 1)(1 ® AC A J + 1 ® 1 ® 1) 
= (1 ® AC A J(1 ® AC A J + (1 ® AC A J(AC An ® 1) 

+(1 ® AC A J(AC An ig) 1)(1 ® AC A J + (AC Aii ® 1)(1 ® AC A J 

+2(1 ® AC A J + (AC An ® 1) + 1 ® 1 ® 1. 

Hence 

(H) 12 (H) 23 (H) 12 - (H) 23 (H) 12 (H) 23 = 1 + 11 + III, (61) 

where 

/ = (AC A „®1)(AC A „®1)-(1® AC A J(1® AC A J, 
// = (AC A „®l)(l(8)AC An )(AC An ®l)-(l(8)AC A J(AC An (8)l)(l® AC A J, 
/// = AC An ® 1-1® AC An . 
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Using ([52p we have 

I = (AC A „®1)(AC A „®1)-(1® AC A J(1® AC A J 
= (AC An ) 2 ® 1 - 1 ® (AC A J 2 

= -(2AC An - n{n + 2)1 ® 1) ® 1 + 1 ® (2AC An - n(n + 2)1 ® 1) 
= -2(AC An <g)l-l<g)AC An ). 

Therefore 

J + // = 1 ® AC An - AC An ® 1. (62) 

By Lemma 2 we get 

J// = (AC An ®l)(l®AC A J(AC An ®l)-(l® AC A J(AC A „®1)(1® AC A J 
= AC An ® 1-1® AC Ari . 
Therefore 

(W)l 2 (W)23(W)l2 - (H)23(H)l2(H)23 =1 + 11 + 111 = 0. 



3.3 Temperley-Lieb Algebraic Structures and Representations 

An L-state TL algebra is described by the elements e*, z = 1, 2, ...,L, satisfying the TL 



algebraic relations [16 



(63) 

CjCj Cj6j , if [z j I ^ 2 , 

and 

= (64) 
where /3 is a complex constant and i = 1, 2, • ■ ■ , L. 

In this section we indicate that there is a TL algebraic structure related to the in- 
tegrable chain model fl59|) , in the sense that the model gives a representation of the TL 
algebra. We suppose that the representation of an L-state TL algebra on an L + 1 chain 
is of the following form, 

d = li ® 1 2 ® • ■ • ® li-i ® £ ® l i+2 ® • • • ® , (65) 
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where 1 is the (n + 1) x (n+1) identity matrix as in section 3.2 and E is a (n + 1) 2 x (n + 1) 2 
matrix. According to formulae (|64l) and (|63"D E should satisfy 



E 2 = (3E . (66) 
(E ® 1)(1 ® £)(£ ® 1) = £ ® 1, 

(67) 

(1® £)(£® 1)(1®£) = 1® £. 



[Theorem 4]. For a given representation of the TL algebra of the form flog) with 
satisfying (|66|) and ( |B"7l) we have that 

fl=E+ -0±vw-* 101 m 



is a solution of the QYBE Q37|). 

[Proof]. For simplicity we set c = (— (3 ± \f(3 2 — 4)/2. Substituting (|68|) into equation 
( |37D and using relations (^) and (p7| ) we get 

R12R23R12 — R23R12R23 

= (E ® 1)(1 ® ® 1) - (1 ® E){E ® 1)(1 ® £7) 

+ c(£ 2 ® 1 - 1 ® E 2 ) + c 2 (E ® 1 - 1 ® E) 
= (£ ® 1)(1 ® £)(£ ® 1) - (1 ® £7)(E ® 1)(1 ® £) + (c/3 + c 2 )(E ® 1 - 1 ® £?) 
= (c/3 + c 2 + 1)(£ ® 1 - 1 ® E) = 0. 

■ 

In general however the converse does not hold, i.e., for a given solution R of the QYBE 
(|3~TD, there does not necessarily exit a TL algebraic representation of the form ( |65D with 
E = aR + b satisfying ( |5"6D and (^) for any constants a and 6. Nevertheless the solutions 
7i of the QYBE in our A n symmetric integrable model fl5T3| ) do give rise to TL algebraic 
representations in the following sense: 

[Theorem 5]. The following {n + l) 2 x (n + l) 2 matrix 

E = + 1 ® 1 (69) 

n+1 



gives the L-state TL algebraic representation (|65D with (3 = 2. 

[Proof]. What we should check is that E in ( |69"P satisfies equations (^) and (j6T|) ■ By 
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Lemma 1 we have 

E 2 = -+I0l 2 = A " \ + 1 1 2 

v n+ 1 7 v n+1 

(AC A J 2 -2nAC An +n 2 l0l 
(n + l) 2 

(-2(n + l)AC An + 2n(n + 1)1 1 
(n + l) 2 

= f3E = 2E, 

i.e., /? = 2. 

From Lemma 1 and fl54f) in Lemma 2 we get 
(#0 ®E)(E® 1) 

= (1®1®1 -)(1 0101 -)(1 0101 - 

n+1 n+1 7i + l 

1 :[(AC An 0l)(l0AC A J(AC An 0l) 



(n + 1) 3 

- n((AC in 1)(1 AC An ) + (10 AC A J(AC An 1)) 

- ti(AC a J 2 1 + 2n 2 AC A „ 1 + n 2 l AC A „ - n 3 l 1 1] 

-^—[(AC An l)(l AC ^J) {AC A 1) 

- n((AC An 1)(1 AC An ) + (10 AC A J(AC An 1)) 

+ 2n(n + l)AC An 1 + n 2 l AC An - 2(n 3 + n 2 )l 1 1] 

, ~ ,,J (n + l) 2 AC An 1 - n(n + 1) 2 1 1 1] 
(n + 1) A 

nl 1 1 - AC A 1 _ „ 

r—r\ — — = ^01. 

(n + 1) 
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By using Lemma 1 and formula Q55D in Lemma 2 we then conclude that 
(1 <g> E)(E (g> 1)(1 ® E) 

= (1 ® 1 ® 1 1 ® 1 <g> 1 -)(1 ® 1 <g> 1 -) 

n + 1 n+1 n+l 

' ■[(l®AC A J(AC An ®l)(l®AC A J 



(n + l) 3 ' 

- n((AC in ® 1)(1 ® AC An ) + (1 ® ACaJ(ACa„ ® 1)) 

- nl ® (AC A J 2 + 2n 2 l <g> AC An + n 2 AC in <g> 1 - n 3 l ® 1 <g> 1] 

' : [(A^ n ®l)(l®AC A J(A^ n ®l) 



(n + l) 3 ' 

- n((AC An ® 1)(1 ® AC A J + (1 ® ACUJ(ACU„ ® 1)) 

+ 2n(n + 1)1 <g> AC An + n 2 AC An ® 1 - 2(n 3 + n 2 )l <g> 1 ® 1] 

— — L-[(n + 1) 2 1 <g> AC A - n(n + 1) 2 1 ® 1 ® 1] 
(ra + l) 3 

nl(g)l®l-l(g) ACa 



1<8>£. 



From fl69|) we see that the Hamiltonian of the A n symmetric integrable chain model 



(p9 ) can be expressed by the TL algebraic elements 

L L L 

H An = EC^W = E( n + + = E( n + ^ + ( n + ( 7 °) 

1 = 1 4 = 1 4 = 1 

with e 4 as in (|65| ) and -E 1 as in (^). Hence instead of the algebraic Bethe Ansatz method, 
the energy spectrum of H An can also be studied by using the properties of the TL algebra 
[ p] ] (for the case of Heisenberg spin chain model, n = 1, see [|TS[]). 

4 Integrable Models and Stationary Markov Chains 
4.1 Stationary Markov Chains 

We first briefly recall some concepts of the theory of Markov chains (for a detailed math- 
ematical description of Markov chains, we refer to [|19|| ). Let Q denote the sample space 
(the set of all possible outcomes of an experiment). If Q is a finite or countably infinite 
sample space and if P is a probability measure defined on the cr-algebra of all subsets of 
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Q, then the pair (Q, P) is called a probability space. A subset A of Q is then called an 
event with probability P(A). 

A function X = X(u), uj G f2, that maps a sample space into the real numbers is 
called a random variable. A stochastic process is a family (X t )tei, I a certain index set, 
of random variables defined on some sample space Q. If I is countable, i.e., I G IV, the 
process is denoted by X 1; X 2 , ... and called a discrete-time process. If I = M + , then the 
process is denoted by {X t }t>o and called a continuous-time process. 

The ranges of X (a subset of real numbers) is called the state space. In what follows 
we consider the case where the state space S is countable or finite. In this case the related 
stochastic process is called a (stochastic or random) chain. 

Let (Q, P) be a probability space in above sense and E, F be two subsets of Q. We 
denote by P(E\F) the (conditional) probability of E given that F has occurred. A 
discrete-time stochastic process {Xi}, % = 1,2,... with state space S = IV is said to 
satisfy the Markov property if for every I and all states ii,t2, it is true that 

P[X t = 2;|AVi = ^_i,X/_ 2 = i z _ 2 ,...,Xi = ii] = P[Xi = U\Xi-i = 

i.e., the values of X_2, X% in no way affect the value of Xi, given the value of X_i. 
Such a discrete-time process is called a Markov chain. It is said to be stationary if the 
probability of going from one state to another is independent of the time at which the 
transition is being made. That is, for all states i and j, 

P[X l = j\X t ^ =i] = P[X l+k = j\X l+k - X = i] 

for k = —(I — 1), — (/ — 2), —1, 0, 1, 2, .... In this case we set pij = P[Xi = j\Xi-i = i] 
and call the transition probability for going from state i to j. 

For a discrete time stationary Markov chain {Xi}, i G IV, with a finite state space 
S = {1, 2, 3, m}, there are m 2 transition probabilities {p^}, i,j = 1, 2, m. P = (p^) 
is called the transition matrix corresponding to the discrete-time stationary Markov chain 
{Xi}. The transition matrix P has the following properties: 

771 

Pij > 0, £]pij = l, i,j = 1,2,... ,m. (71) 

i=i 



Any square matrix that satisfies condition (71) is called a stochastic matrix 



24 



A continuous-time stochastic process, {X t }t£M + is said to satisfy the Markov property 
if for all times to < t\ < ... < ti < t and for all I it is true that 

P[X t = j\X to =i ,X tl =i 1 ,...,X tl =i l ]=P[X t = j\X il =i l ]. 

Such a process is called a continuous-time Markov chain. It is said to be stationary if for 
every % and j the transition function, P[X t +h = j\Xt — i], is independent of t. In this 
case P(t) = P(X t= j\X = i) is a semigroup (e.g. on l 2 (S)), called transition semigroup 
associated with the Markov chain. Its generator Q = (g^) has the properties: 

Qij > 0, i^j, qu = -J2 ( liii i,j = l,2,...,m (72) 

and is called an intensity matrix. Vice versa, any Q (satisfying ( [72]) and properly defined 
as a closed operator when S is infinite) gives rise to a unique continuous-time transition 
semigroup, P(t) = e™, t > 0, which can be interpreted as transition semigroup associated 
to a certain Markov chain (with state space S) ||19|| . 

The properties of Markov chains are determined by the transition matrix P for 
discrete-time stochastic process and the intensity matrix Q for continuous-time stochastic 
processes. If the eigenvalues and eigenstates of P and Q are known, then exact results 
related to the stochastic processes, such as time-dependent averages and correlations, can 
be obtained. 

Now we consider a chain (in the algebraic sense of sections 1-3) with L + 1 sites. To 
every site i of the chain we associate n+ 1 states described by the variable Tj taking n+ 1 
integer values, 

n = (7f = 0, r] = 1, Tf = 2, t? = n), (73) 

(conventionally a vacancy at site i is associated with the state 0). We associate to the 
lattice site i of the algebraic chain a Hilbert space of dimension n + 1. The state space of 
the algebraic chain is then finite and has a total of (n + 1) L+1 states. 

By definition, for an integrable chain model with Hamiltonian H one has exact solu- 
tions for the eigenvalues and eigenstates of H. The system remains integrable if one adds 
to H a constant term c and multiplies if by a constant factor d . Moreover the eigenvalues 
of H will not be changed if one changes the local basis, i.e., the following Hamiltonian 
H', 

H' = BHB-\ B = ®£?Bi, (74) 
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where Bi = b and b is an (n + 1) x (n + 1) non singular matrix, has the same eigenvalues 
as H. Therefore if an integrable chain model with Hamiltonian H can be transformed by 
B (modulo c, d) into a stochastic matrix P, in the sense that 

P = B(c'H + cl)B-\ (75) 

where 1 is the (n+ x (n + identity matrix, B as in (fT4[), such that P satisfies 
(|7lD (with m = (n + 1) L+1 ), then P defines a discrete-time Markov chain and the related 
stochastic process can be simply studied by using the properties of the related integrable 
model with Hamiltonian H . 

And if an integrable chain model with Hamiltonian H can be transformed (modulo c, 
c') into an intensity matrix Q, in the sense that 

Q = B(c'H + cH)B- 1 (76) 

with Q satisfying (|72"D (with m = (n + 1) L+1 ), then Q determines a continuous-time 
Markov chain and its properties can also be obtained by using the results of the related 
integrable Hamiltonian H. 

In the following we discuss the question of wether the integrable models obtained 
in the way presented in this paper could be transformed into stationary Markov chains 
through transformations of the forms ( [75]) or (|7(J). Of cause a stochastic matrix P (i.e. as 
given by ([H])) can not in general be transformed into a intensity matrix Q (as given by 
(|72"D) by the spectrum preserved similarity transformation, P = B(d Q + c\)B~ l . That is, 
an integrable chain model with Hamiltonian H that gives rise to a discrete-time Markov 
chain by the transformation ( |75| ) will in general not give rise to a continuous time Markov 
chain by ([76]), and vice versa. 

4.2 Discrete-time Markov Chains Related to A n Symmetric In- 
tegrable Models 

We first note that for an integrable chain model with Hamiltonian H = J2i=i hi.i+i and 
(n + 1) states at every site i, i = 1, 2, L + 1, if the sum of the elements in any row of 
the (n + l) 2 x {n + l) 2 matrix h is 1/L, then the sum of the elements in any row of the 
matrix if is 1. Hence if under the following transformation h — >• h! given by 

h! = {b® b)(c'h + cl <g> 1)(6" 1 ® r 1 ), (77) 
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the sum of the elements in any row of hi is 1/L and (h') a ,p > 0, a, (5 = 1, 2, (n + 1) 2 , for 
some real constants c', c and a non singular (n+1) x (n+1) matrix 6, then P = X)i=i ^4 j+i 
defines a stationary discrete-time Markov chain. P has the same eigenvalue spectrum 
(shifted by a constant) as the spectrum of the integrable model with Hamiltonian H. If 
P is invariant under a certain algebra A, we call the Markov chain A symmetric. 

[Theorem 6]. The following matrix 

L 



1 



L(n + 1) 



H 



1 



L{n+l)U 



£(Ac An + 1 ® i) w 



E 



n(n+l)/2 



a=l 



L{n + 1) ^ 

+ a ( n + 1 - a)(^a)i(/ia)i+l 



a=l 
n n—a 



+ EE a ( n + 1 - " - P)((h a )i(h a+f3 ) i+1 + {h a+/3 )i(h a ) i+1 ) 

a=l 0=1 

defines a stationary discrete-time A n symmetric Markov chain. 
[Proof]. I. Set h' = I ^i rTy (AC An + 1 <g> 1). Then 

L 



+ 



n + 1) 



(7S 



(79) 



From formula 



we have 



{h>) 



1 



a/3 



L(n + 1) 
1 



(Ac An + 1 ® i; 



a/3 



[5 Q ,s[(n + l)5 a ,z( n+ l) + z + l - 1] 



L(n + 1) 

+ (n + l)[5 a j( n+ 2)+fe+2^/3,(i+l)(n+2)+A;(n+l) 
+^/3J(n+2)+fc+2^Q,(j+l)(n+2)+fc(n+l)] + ^a/j] 

= -j\^ a ^ a ^ n+1 )+ l+1 ^aJ(«+2)+fc+2^/3,(j+l)(n+2)+fc(n+l) 
+^/3J(n+2)+fc+2^a,(j'+l)(n+2)+fc(n+l)] > 0. 

Therefore (P An ) a p > 0, a, (3 = 1, 2, (n + l) 2 . 



(80) 
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II. By using the identity (|53|), we get 

(n+l) 2 , (n+1) 2 

X! = 77 — TT H + l)<W<W(n+l)+Z+l 

/3=1 M 71 + i J 0=1 

+ (n + l)[5 Q j( n +2)+fc+2^/3,(j+l)(n+2)+fc(n+l) 
+^/3J(n+2)+fc+2^Q,(i+l)(n+2)+fc(n+l)]] 

1 1 

(n + l) = -. 



L(n + l) v ' L' 

Hence the sum of the elements of any row of the matrix P^ n is one, i.e., Y^b=P (Pa*) at 

v (n+l) i+1 / V i 7/ \ i 

2^/3=1 lZ^j=l ,l 'i,i+l)aP — L - 



III. As is invariant under A n , P^ n = jrj^\ is obviously invariant under A n and 
has the same spectrum as H^ n - 

By the definition (|7lD Pa„ is the transition matrix of a stationary discrete-time A n 
symmetric Markov chain. ■ 

We give some discussions on the properties of the stationary discrete-time A n sym- 
metric Markov chain associated with the stochastic matrix Pa u - The state space of this 
Markov chain is, S = (1, 2, (n + 1) L+1 ), which corresponds to (n + 1) L+1 states, 

(r ®TL® ...®r n ), (81) 

Tj as in (|73|) , of the algebraic chain with L + 1 lattice sites. 

The properties of a Markov chain are determined by the transition matrix P = (pij). 
A subset C of the state space S is called closed if = for all i G C and j ^ C. If a 
closed set consists of a single state, then that state is called an absorbing state. A Markov 
chain is called irreducible if there exists no nonempty closed set other than S itself. 

From formula ( |5(iD we have 

{h')aa = (^')(n+l) 2 ,(n+l) 2 = X> (h')a/3 = (h')/3a = 0, (3 ^ a, 

a = l(n + l) +1 + 1, / = 0,l,...,n. 

Let 



(82) 



o ( I > + l)((n + l) L -l)+n \ 
oo = a a = t HI, / = U, 1, n, (8d) 

be a subset of the state space S. From formula (|79|) , with 

(/i')m+i = li ® I2 ® ••• ® l<-i ® ft' ® li+a ® ••• ® 
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we get 

(P An ) aa = l, (P An )(3 a = (PA n ) al 3 = 0, P^a, aeS . (84) 

Therefore the n + 1 states in Sq are absorbing states of the Markov chain Pa„- This 
chain is by definition reducible. For a reducible Markov chain the "long time" probability 
distribution, if it exists, may depend on the initial conditions, i.e., lim^ 00 (P J 4 n )^ may 
depend on 7. From the properties ( ERj ) of Pa„, we see that if the Markov chain P An 
is initially at one of the states a £ Sq, it will remain in that state a forever. These 
n + 1 absorbing states correspond to the states of the algebraic chain through (|8~i"|). For 
instance, the states 1 and (n+ 1) L+1 in S correspond to the states (0, 0, 0) (all the sites 
of the algebraic chain are at state 0) and (n, n, ...,n) (all the sites of the algebraic chain 
are at state n). 

4.3 Continuous- time Markov Chains Related to A n Symmetric 
Integrable Models 

For an integrable chain model with Hamiltonian H = J2f=i arid with (n + 1) states 
at every site of the chain, if the sum of the elements in any column of the matrix h is 0, 
then the sum of the elements in any column of the matrix H is also 0. Hence if under the 
following transformation h — > h" with: 

h" = (b ® b)(c'h + cl® ® 6 -1 ), (85) 

the sum of the elements in any column of h" is and (h") a ^ > 0, a ^ (3 = 1, 2, (n+ 1) 2 , 
for some real constants c', c and a non singular (n + 1) x (n + 1) matrix b, then Q = 
J2i=i h'ii+i i s the intensity matrix for some stationary continuous-time Markov chain. 
Q has the same eigenvalue spectrum (shifted by a constant) as the spectrum of the 
Hamiltonian H. We call the Markov chain A symmetric if Q is invariant under the 
algebra A. 
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[Theorem 7]. The following matrix Q is the intensity matrix of a stationary continuous- 
time Markov chain, 



L 

Q An = H An -{n + l)L = Y J (^C An -nl®l) hi+l 

i=l 

n(n+l)/2 n 

(n+1) J2 i( e a)i(f a )i+i + (f a )i(e Q )i+i) + a ( n + 1 ~~ a ){K)i{K)i+l 



E 

1=1 



a=l a=l 
n n—a 



+ EE a ( n + 1 - " - P)((h a )i(h a+/3 ) i+1 + (h a+ p)i(h a ) i+1 ) 

a=l (3=1 



— nL. 

(86) 



[Proof]. Set h" = (AC An - nl ® 1). Then 

<5a„ =Y / h'li+v (87) 



8=1 



From ( [48|) we observe that, for a ^ (3, 

h'a-t/3 = {n+ l)[5aj(ri+2)+fc+2^/3,(j+l)(n+2)+fc(n+l) 
+^/3j'(™+2)+fc+2^a,0'+l)(n+2)+fc(n+l)] > 0. 

Therefore (Qa„ W/3 > 0, a, /3 = 1, 2, .., (n + 1) L+1 . 

Again by (|48|) the sum of the elements in any given column (5 of the matrix h" is 



(n+l) 2 (n+1) 2 

E h 'Lp= ( n + 1 )[<M ( W(n+l)+J+l - 1) 

a=l a=l 

+ (n + l)[5 a ,j(n+2)+A:+25/3,(j+l)(n+2)+fc(7i+l) + ^J(n+2)+fc+2#a,(j+l)(n+2)+fc(n+l)] 
(n+1) 2 

= ^2 (n + l)[—8 a p + 5 a ,j(n+2)+A:+2^/3,(j+l)(n+2)+fe(n+l) 

a^(n+l)+Z+l 

+ ^/3,j(n+2)+JH2^a,0'+l)(n+2)+ft(n+l)j 
= (n + 1)( — 1 + ^/3,(i+l)(n+2)+fc(n+l)U=j(n+2)+fc+2 + <5/3,j(n+2)+fc+2 1 a=(j+l)(n+2)+fc(n+l) ) 

= o, 

i.e., the sum of the elements in any given column /3, (3 = 1, 2, (n+ l) 2 , of the matrix /i" 
is zero. Therefore the sum of the elements in any given column /3, (3 = 1, 2, (n+ 1) L+1 , 



of the matrix Q An is also zero, Ea=i " (QA„)a/3 = 0. At last Q An = H An - (n + l)L i 



is 



obviously A n symmetric with the same spectrum (shifted by a constant) as H An . ■ 

The long run distribution of the Markov chain described in Theorem 7 is given by the 
vector 7r = (tti, 7r 2 , ...), where 7Tj represents the "long time" probability of the state i G S, 
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satisfying 

0+l) L+1 (n+l) L+1 

E (QA n )apir a = 0, V/3gS, J2 7r « = 1 - ( 88 ) 

However as this Markov chain is not irreducible, the solution of the equation ([88]) is not 
unique but depends on the initial conditions. From (|48|) we see that 

(h")ap = (h") Pa = 0, V/3, 

a = l(n + l) + l + l, / = 0,l,...,n. 

Hence from (|87|) we get 

(QaJoP = (QA n )p a = 0, a G So, V/3, 

with So as in (^). Therefore if this Markov chain is initially at a given state a G So, it 
will remain at that state. 

The states (3 $ S form a closed subset of S. From (|48| ) and (|8T| ) one also learns 
that the absolute value of all the nonzero elements of any column of the intensity matrix 
Qa„ are equal. Let S' be a closed subset of S with I elements. If the Markov chain is 
initially in the closed set S', then it will remain in S' and the long run distribution is 
7r = (jn, 7r 2 , 7T( n+ i)L+i), where 7Tj = l/l for i G S' and 7Tj = if i £ S' . 

5 Conclusion and Remark 

Using the Casimir operators and coproduct operations of algebras, we have given a sim- 
ple way to construct chain models with a certain algebraic symmetry and nearest or 
non-nearest neighbours interactions. We discussed integrable chain models with nearest 
neighbours interactions with the symmetries provided by the fundamental representation 
of the classical Lie algebra A n . It is shown that corresponding to these A n symmetric inte- 
grable chain models there are exactly solvable stationary discrete-time (resp. continuous- 
time) Markov chains whose spectra of the transition matrices (resp. intensity matrices) 
are the same as the ones of the corresponding integrable models. 

Other symmetric integrable models ( e.g. with B n , C n , D n symmetry) and related TL 
algebraic structures and Markov chains can be investigated in a similar way. The discus- 
sion of integrable models related to higher dimensional representations of the algebras and 
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the integrability of chain models with non-nearest neighbours interactions is postponed 
to further work. 
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